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Introduction 


Control  of  the  supercavity  parameters  is  necessary  for  maintenance  of  calculated  regime  of  the 
supercavitation  flow  in  unsteady  conditions:  at  changing  the  velocity  and  depth  of  the  body  motion.  The 
problem  of  cavitation  flow  control  includes  the  two  problems  connected  among  themselves: 

I)  the  control  of  cavity  dimensions; 

II)  the  control  of  forces  on  a cavitating  body. 

We  consider  sequentially  both  problems. 

I)  It  is  known  [1,2]  that  at  small  G the  mid-section  diameter  Dc  and  length  Lc  of  the  axisymmetrical 
cavity  are  determined  by  the  cavitator  diameter  Dn  , the  drag  coefficient  Cx  and  the  cavitation 
number: 
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where  the  cavitation  number  <7  = 2 {p  g h — Pc)/  pv2  is  a an  important  parameter  of  supercavitating 
flows.  Here,  h,  V are  the  depth  and  the  speed  of  motion,  respectively,  and  pc  is  the  cavity  pressure. 

At  small  cavitation  numbers  [6]  the  following  approximate  formula  for  the  drag  coefficient  of  a 
blunted  cavitator  is  valid: 


cx=cx0(1  + g),  (2) 

where  cr0  is  the  drag  coefficient  value  for  a given  cavitator  when  <7  = 0 (that  is,  for  free  streamline 
flow). 

The  validity  of  asymptotic  relations  (1)  for  models  with  disk  cavitator  moving  with  velocities 
300  = 1300  m/s  is  confirmed  by  us  experimentally  [3]. 

It  is  visible  from  relations  (1)  that  at  constant  Dn  ,h,V  to  control  cavity  dimensions  it  is  possible  by 
two  ways: 

1)  changing  the  cavitation  number  <7  ; 2)  changing  the  cavitator  drag  coefficient  cx . 


Paper  presented  at  the  RTOAVT  Lecture  Series  on  "Supercavitating  Flows  ”,  held  at  the  von  Karmctn 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 
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1.  Supercavitation  flow  control  by  gas  blowing 


The  first  way  is  effective  at  motion  velocities  10  4- 1 00  m/s  and  is  realized  by  gas  supply  to  the  cavity, 
i.e.  by  increase  pc  (so  named  the  artificial  cavitation  or  ventilation  [5]).  The  method  of  gas-supply  is 
widely  applied  in  experimental  researches  for  formation  of  long  cavities  at  moderate  flow  speeds  in  the 
hydrodynamic  tunnel  (Fig. l).The  Fig.  2 demonstrate  dependencies  0 = g(q  ). 

The  value  of  gas-supply  to  the  cavity  Q is  defined  by  dependence  [2,  6]: 

Q=-^j  = F(0,Fr,Re,We),  (3) 

VDl 

where  Fr,  Re,  We  are  accordingly  Fronde,  Reynolds  and  Weber  numbers.  The  Fig.  3 demonstrate  the 
cavity  length  - L control  by  gas  flow  rate  Q . In  the  region  of  strong  influence  of  gravity  the  formula  for 
supply  value  has  the  form  [2]: 


0 = 


0.27 

a[a3Fr4 -2] 


(4) 


In  the  regimes  close  to  vapor  ones  the  structure  of  dependence  (4)  changes: 


Q=kVS\ 
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(5) 


The  region  of  effective  control  of  air-supply  in  the  artificial  cavitation  regime  is  limited  by  value 
gDc!  V 2 from  (4)  and  by  vapor  cavitation  occurrence  at  0 = 0V : 


1 .4 \gDc  / V2  < a < 0V  = 2(P  Jv\  (6) 

pv 

where  pv  is  the  pressure  of  saturated  water  vapor. 

The  motion  in  water  with  velocities  V »100  m/s  occurs  in  the  regime  of  vapor  or  natural 
cavitation,  then  the  gas-supply  to  the  cavity  becomes  inefficient. 


2.  Supercavitation  flow  control  by  jet  cavitator 


According  to  the  hydrodynamic  scheme  of  supercavitation  flow,  the  object  is  placed  partially  or  fully 
inside  a supercavity  (Fig.  4a,  b)  formed  by  the  nasal  part  (cavitator)  [1,  4,  6].  In  the  case  of  a jet  cavitator 
system,  the  cavity  separates  from  the  craft  hull  (Fig.  4c). 
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The  comparatively  new  scheme  (c)  has  in  comparison  with  conventional  schemes  (a)  and  (b)  the 
following  advantages: 

1)  the  hydrodynamic  drag  reduction  due  to  the  model  in  general  has  not  points  of  contact  with  water 
(theoretically  in  twice  in  comparison  with  a solid  disk  cavitator  at  0 — > 0); 

2)  the  possibility  of  control  of  cavity  by  water  rate  change  in  the  nose  jet. 

For  the  first  time  L.l. Sedov  offered  the  way  of  formation  of  cavities  due  to  giving  up  the  fluid  jet 

toward  to  the  mainstream  ( i.e.  scheme  of  jet  cavitator)  [4],  He  gave  the  estimation  of  drag  coefficient 
of  the  jet  cavitator,  considering  it  as  an  inverted  scheme  of  cavity  with  reentrant  jet  [14], 

The  scheme  of  cavitating  flow  around  the  body  A with  reentrant  jet  is  presented  in  Fig.  5a.  It  is 
known  as  the  Efros  - Hilbarg  scheme  [14],  It  is  theoretically  possible  to  think  that  the  reentrant  jet  is 
absorbed  by  fictitious  body  B . Then  the  cavitation  drag  of  body  A is  equal  to: 


x = pQ(v„  + vc),  vc  = v„y[I7^,  (7) 

where  Q = SVC  is  the  volumetric  water  rate  in  the  reentrant  jet,  S is  the  area  of  jet  section,  Vc  is  the 
velocity  on  the  cavity  surface  (it  is  equal  to  velocity  in  the  reentrant  jet). 

Passing  on  to  the  limit  at  o — > 0 (thus  the  body  B tends  to  the  infinity ),  we  obtain  the  cavitating 
drag  of  body  A by  Kir  chh  off  scheme: 


X0=2pQVoo. 


(8) 


Changing  the  direction  of  all  the  velocities  at  <7  — > 0,  we  obtain  the  scheme  of  cavitating  flow  around 
the  body  B with  jet  cavitator  ( Fig.  5b).  The  expenditure  of  energy  expending  by  body  on  overcoming  of 
cavitator  drag  is  replaced  by  expenditure  of  energy  on  formation  of  opposing  jet  in  this  case.  We  obtain 
for  jet  momentum: 


PQV~=^Y’  (9) 

i.e.  when  0 — > 0 , the  opposing  jet  momentum  necessary  for  cavity  supporting  is  twice  less  than 
cavitating  drag  force  of  a solid  cavitator. 

We  experimentally  investigated  the  cavities  formed  by  opposing  liquid  jet  in  the  small  hydrodynamic 
tunnel  at  the  1HM  UNAS  at  flow  velocity  10  m/s.  Thus  the  air  is  in  addition  supplied  to  the  cavity 
because  of  small  flow  velocity.  The  experimental  dependence  of  opposing  jet  length  from  the 

nozzle  to  the  frontal  point  of  cavity  on  velocity  in  the  jet  V - = V - /¥„  is  presented  in  Fig.  6. 

We  note  that  the  limit  value  obtained  experimentally  Vj  (Lj)  ~ 0.7 at  L , — > 0 does  not 

2 2 

correspond  to  the  classical  scheme  of  jet  collision  with  critical  point  [14],  in  which  pV M = pV j . The 
value  Vj  (0)  ~ 0.7 SV^  permits  to  assume  that  the  scheme  of  jet  collision  with  formation  of  critical 
region  without  breakdown  points  [14]  is  realized  in  the  experiment.  For  such  scheme  the  relation 
pV^=2pVj  should  be  fulfilled  theoretically,  whence  Vj  = 0.7  01 . 
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We  mention  yet  for  completeness  about  a cavitating  flow  scheme  similar  externally  to  the  scheme 
where  the  cavity  is  formed  with  help  of  gas  jet  blown  out  from  the  model  nose  toward  mainstream  [5]. 
Thus  these  cavities  relate  to  artificial  (ventilated)  cavities  by  classifications.  To  create  cavity  the  gas  jet 
with  density  should  be  blown  out  with  velocity  V enabling  to  overcome  the  dynamic  head  of 

water  mainstream  pV^  12  . We  have  from  the  Bernoulli  integral: 


Vg>Voa  -?-{\  + o). 
V Pg 


(10) 


If  the  air  is  used  as  a gas  p / pa  = 800 . then  we  obtain  from  (10)  without  account  of  air  compressibility: 


V„  >28.3  1 + — 


K . 


) 


It  is  visible  that  such  way  of  cavity  formation  is  really  realized  only  at  small  Vx  . 

The  conducted  experiments  show  that  the  cavities  created  due  to  opposing  gas  jet  have  the  strongly 
disturbed  surface  and  the  powerful  accompanying  turbulent  emulsion  flow  [5], 


3.  Supercavitation  flow  control  by  change  of  cx 


As  it  is  was  said  above  the  cavity  control  method  by  means  of  gas-supply  is  not  applied  at  very  high 
velocities  of  motion  in  water  and/or  small  model  dimensions.  Therefore,  the  research  of  possibility  of 
cavity  control  by  change  of  C is  very  interesting. 


3.1.  Estimation  of  drag  of  axisymmetric  cavitator 


The  practical  way  of  estimation  of  the  drag  of  the  axisymmetric  cavitators  consists  in  that  the 
pressure  distribution  p along  the  generator  of  cavitator  is  taken  from  the  solution  of  the  two- 
dimensional  problem  on  free  streamline  flow  around  the  contour  coinciding  with  a meridian  section  of 
the  cavitator.  Then  the  drag  of  the  axisymmetric  cavitator  is  defined  by  integration: 


c,  = 


2X 

pv2kR2„ 


X = 2k  j(p- pjydy. 


(11) 


where  R„  = Z)  / 2 . 
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It  was  shown  in  the  work  [17]  that  the  values  Cx  for  cones  computed  by  given  approximate  way  are 

close  to  the  experimental  data.  We  present  also  their  comparisons  with  results  of  the  numerical 
computation  executed  on  basis  of  exact  formulation  of  the  problem  [18].  For  comparison  we  use  the 
approximation  formulae  given  in  the  work  [19]: 


cx  = cx o + (0.524  + 0.672Jtf)<r,0  < a < 0.25, 

1 1 

12  2 

cx0  = 0.5  + 1.81(|/  - 0.25)  - 2{ji  - 0.25)2, 

1 1 

-<H<- 
12  2 


cx 0 =^(0.915  + 9.5^),  0</i<—, 


(12) 


where  /un  is  the  half  angle  of  cone.  When  <7  — > 0 , the  values  Cx0  are  obtained  by  extrapolation 

method.  A comparison  of  the  experimental  data  with  the  results  of  "exact"  [19]  and  approximate  [17] 
computations  is  given  in  Table  1 : 


Table  1 


rtn  (deg.) 

15 

30 

45 

60 

75 

90 

Exper. 

0.15 

0.35 

0.47 

0.61 

0.72 

0.82 

[7] 

0.1428 

0.3353 

0.5000 

0.6369 

0.7461 

0.8275 

[io] 

0.2045 

0.3758 

0.5181 

0.6350 

0.7296 

0.8053 

We  obtain  the  calculated  formula  for  an  arbitrary  contour  from  (11): 


cx0 


(13) 
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3.2.  Cavitator  with  variable  drag 


Both  the  length  and  the  biggest  diameter  of  cavity  are  proportional  to  ^cxo  as  it  is  seen  from  the 
relations  ( 1 ). 


D^_  ko  (1+<?) 

Dn  V cr 


A 


1 


D„ 


<J 


1 


^ = - cx0(l  + a)\n- 


(7 


(14) 


Hence,  the  cavity  shape  should  change  similarly  to  self  at  changing  the  control  parameter  cx 0 when 

Dn  and  (T  are  constant.  It  may  be  obtained  by  changing  the  cavitator  shape  at  constant  diameter  of  the 
cavitating  edge. 

We  offer  the  sketch  of  the  axisymmetric  cavitator  with  variable  geometry.  The  meridian  section  of 
the  cavitator  has  the  £ - shaped  contour.  The  cavity  drag  coefficient  is  increased  from  cx  ( /3)  up  to 
1 for  a cone  at  moving  of  the  internal  cone  relatively  to  the  outside  cartridge  with  a cavitating  edge,  i.e. 
at  increase  of  parameter  x / Dn  from  0 . 

The  graphic  of  dependence  cx(i  on  x / Dn  is  shown  in  Fig.  7 for  ^ - shaped  contours. 

The  graphic  of  dependence  of  the  ratio  cv0  for  axisymmetrical  cavitators  to  theoretical  value  of  the 

drag  coefficient  for  disk  0.8053  on  the  dimensionless  extension  of  cavitating  edge  x/ Dn  is  shown  in 

Fig.  8.  Also,  the  experimental  points  are  plotted  there.  As  it  is  seen,  the  drag  increases  slower  in  the 
experiment  than  as  the  computation  predicts  at  increase  of  x . It  is  explained  that  the  closed  zones  of 
vortical  motion  of  water  are  formed  in  the  interior  angles  of  cavitator  in  the  experiment.  Other  words,  the 
potential  scheme  of  the  free  streamline  flow  with  the  critical  point  is  not  realized  in  this  case,  but  the 
scheme  of  the  flow  with  stagnant  zones  is  realized  [14]. 

The  appearance  of  cavity  at  the  two  extreme  positions  of  cone  (d  = 0.071,(1  = 1/3)  is  given  in  Fig. 
9.  At  extending  of  cavitating  edge  from  x / Dn  = 0 to  x/  Dn  =0.138  in  the  experiment  we  obtained  the 
increase  of  main  dimensions  of  cavity  on  35  %. 


4.  Stability  of  supercavitating  motion  of  bodies 


4.1.  Possible  flow  schemes 

The  analysis  showed  that  the  four  different  mechanisms  of  motion  stabilization  sequentially  act  at 
motion  velocity  increase. 

1.  Two-cavity  flow  scheme  (Fig.  10a),  V ~ 70  m/s. 

In  this  case  the  hydrodynamic  drag  center  is  placed  behind  the  mass  center,  and  stabilizing  moment  of 
the  force  Y2  acts  to  the  model.  It  means  that  the  classic  condition  of  stability  is  fulfilled. 
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2.  Stationary  planing  along  the  internal  surface  of  cavity  (Fig.  10b),  V~  50  -4-200  m/s.  To 
compensate  the  buoyancy  losses  the  body's  tail  part  is  planing  along  the  lower  internal  cavity  surface. 
In  this  case  the  motion  may  be  stable  as  a whole. 

3.  Impact  interaction  with  cavity  boundaries  (Fig.  10c),  V ~ 300  -4-900  m/s.  Presence  of  initial 
perturbations  of  the  model  attack  angle  and  the  angular  velocity  causes  to  the  impact  of  the  model  tail 
part  against  the  internal  boundary  of  cavity.  The  mathematical  modeling  showed  that  after  this  impact 
the  model  can  perform  steady  or  damped  oscillations  accompanied  by  periodic  impacts  of  the  tail 
part  alternately  against  the  upper  and  lower  cavity  walls.  Then  the  motion  can  remain  stable  as  a whole. 

4.  Aerodynamic  interaction  with  vapour-splash  medium  of  cavity  (Fig.  lOd),  V ~1000  m/s  and 
higher.  At  very  high  velocities  the  aerodynamic  and  splash  forces  of  interaction  with  vapour  filling 
the  cavity  and  splashes  near  the  cavity  boundaries  have  considerable  effect  on  the  body  motion. 
Since  the  clearance  between  the  body  surface  and  the  cavity  boundary  usually  is  small  compared  to  the 
cavity  radius,  we  apply  the  known  methods  of  near-wall  aerodynamics  to  estimate  the  arising 
forces. 


Statically  stable  contours 

The  pressure  forces  acting  on  the  inclined  contours  in  free  streamline  flow  sum  to  the  resultant  force 
F 

and  moment  M0  showninFig.il. 

The  contour  will  be  statically  stable  if  these  forces  tend  to  turn  the  cavitator  to  decrease  angle  a 
relative  to  the  mass  center  of  the  model.  The  condition  of  static  stability  has  the  form: 

Fv 

<5=arctan—  >a,  (15) 

F 

1 X 


where  Fx  and  Fv  are  the  projections  of  the  vector  F onto  the  x-y-axes. 

Computations  have  shown  that  a cavitator  is  statically  stable  if  its  contour  is  concave  towards  to  the 
stream  ( in  particular,  inverted  wedges  with  [l  > 0.5 ).  A flat  plate  will  be  neutrally  stable, 

Fy  / Fx  = tana , if  we  do  not  take  into  account  the  moment  Mq  induced  by  displacement  of  point  of 

application  of  force  F from  the  point  C.  Computations  have  shown  that  this  moment  is  always 
stabilizing,  but  is  very  small  in  value.  Convex  contours  ( in  particular,  wedges  with  at  jd  < 0.5  ) are 
statically  unstable. 

The  force  polars  for  wedges  are  shown  in  Fig.  12  for  a series  of  values  of  wedge  semiangle  {3  =Kjd 
in  degrees.  The  characteristic  stability  regions  are  shown  in  Fig.  13.  Wedges  are  statically  stable 
for  / 3 > 90° , and  are  statically  unstable  for  [3  < 90“ . The  lift  of  a wedge  is  equal  to  0 irrespective  of 
the  attack  angle  for  [3  ~ 50°35'  [9]. 

Computer  program  FLOWJET  allows  to  select  the  best  cavitator  shape  for  static  stability  when  the 
arbitrary  is  maximum.  The  best  cavitator  from  a practical  point  of  view:  (1)  is  statically  stable;  (2) 

has  a maximal  rotational  derivative  ; and  (3)  has  a minimal  drag  coefficient  cx0  . 
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4.2.  Designing  the  supercavitating  (SC)  models 


The  supercavitation  flow  scheme  is  a special  artificial  scheme.  It  is  necessary  to  know  a cavity  shape 
for  successful  realization  of  this  scheme.  This  is  necessary  for  optimal  location  of  the  body  in  the 
supercavity. 

We  consider  two  problems: 

- calculation  of  a supercavity  for  a body  of  given  shape  and  given  motion  regime; 

- design  of  optimal  body  shape  for  the  known  supercavity. 

Stage  I.  Inscribing  the  body  into  the  cavitv 

According  to  the  given  motion  regime  - V = V(t);  H = H(t)  we  calculate  minimal  and  maximal 
cavitation  numbers 


^min  ’ ^min  ^ ^max  ^K^min  > ^max  ) 


to  construct  minimal  and  maximal  supercavity  contours  R = F(Rn;  cf)  past  a circular  disk  [16]. 

The  main  purpose  of  this  stage  is  to  inscribe  the  model  into  the  minimal  cavity  contour.  The  minimal 
contour  is  determined  by  presence  of  points  of  contact  of  the  model  contour  and  the  free  cavity 
boundary  in  the  tail  part. 

Fig.  14  presents  results  of  the  computation  of  cavity  shape  past  disk  cavitators  around  the  body  of 
revolution.  These  results  are  obtained  with  the  help  of  software  SCAV  developed  at  IHM  UNAS.  The 
program  permits  to  plot  supercavities  by  given  parameters:  diameter  and  angle  of  the  cavitator  Dn ; [3  ; 
velocity  V and  depth  II  of  the  motion:  the  pressure  in  the  cavity  P,  . 

At  that  we  determine  the  clearance  h between  the  body  surface  and  the  cavity  boundary,  ring  area  of 
this  clearance  along  the  whole  body  length,  moment  of  inertia  and  body  mass  (Fig.  15). 

The  program  also  computes  the  distance  and  the  velocity  reduction  from  the  starting  velocity  V = V(l 
upon  motion  on  inertia  and  displays  graphs  of  changing  the  length  Lc  and  diameter  D of  the  supercavity 
as  well. 


Stage  II.  Optimization  of  the  supercavitating  body  contour 

This  stage  purpose  is  to  optimize  the  clearance  between  the  free  cavity  boundary  and  model. 

During  motion  of  the  model  by  given  regime  (V  = V(l\  II  = //(/))  the  supercavity  shape  will 
change  according  to  changing  the  cavitation  number 


^ min  <0<0 


max  • 


The  clearance  between  the  body  and  cavity  boundary  will  be  changed  in  certain  range  according  to 
changing  of  <7. 
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Besides  the  clearance  will  be  changed  under  action  of  possible  angular  oscillations  of  the  model  about 
the  center  of  mass  in  supercavities  with  dimensions  exceeding  the  minimal  possible  ones  (c  > crmin ). 

Fig.  16  shows  a PC-screen  copy  of  optimization  procedure  for  the  model  shape  at  its  angular 
inclination  y/.  Simultaneously,  values  of  upper  h2  and  lower  h\  clearances  between  the  model  and 
supercavity  contours  are  shown  there.  The  program  STAB  permits  to  determine  the  maximal  angles  of 
trim  y/  = l//max , points  of  touching  the  model  with  the  cavity  and  the  mass  center  location.  A simple  rule 
was  developed  for  ensuring  the  stability  of  the  supercavitating  model  motion  at  the  1HM  UNAS. 


Rule  of  optimization  of  SC  model  shape 


For  certain  ranges  of  varying  the  cavitation  number  dmin  < <7  < CTmax  and  angles  of  trim  i//  < |i//max| , 

the  optimal  shape  of  the  supercavitating  model  must  not  have  points  of  touching  with  the  supercavity 
contour  which  are  situated  in  front  of  the  model  center  of  mass. 

This  requirement  does  not  take  into  account  the  dynamics  of  the  body  and  is  based  on  principles  of 
static  stability.  In  the  case  of  continuous  flow  the  analogous  rule  of  the  stable  motion  maintenance 
requires  location  of  a point  of  the  hydrodynamic  force  application  past  the  body  center  of  mass.  Both  the 
rules  guarantee  against  appearance  of  destabilizing  moments  directed  to  the  possible  angle  of  model 
deflection.  This  means  that  the  induced  total  moment  acting  on  the  model  always  must  has  a sign 
opposite  to  a sign  of  angle  y/  - accidental  deflection  of  the  model 

SgnM  = -Sgn !//  . 


4.3.  Equations  of  the  SCM  dynamics 

The  most  effective  method  of  investigation  of  the  supercavitating  body  dynamics  proved  to  be  its 
direct  simulation  on  a PC-screen  (computer  experiment).  It  permits  to  research  the  motion  stability  in 
interactive  regime  "researcher  - computer"  [16]. 

The  complete  mathematical  model  of  the  SCM  motion  includes  a set  of  equations  of  solid  body 
dynamics,  equations  to  calculate  the  unsteady  cavity  shape  and  relations  to  calculate  the  acting  forces. 

A set  of  dynamic  equations  of  the  axisymmetric  body  motion  in  vertical  plane  without  rotation  in  the 
body  co-ordinate  system  OjXjjq  is 


m\ 


dV 

dt 


+ cdxF 


ZF, 


(16) 


) 


T d ® vCf 

L = ZM, 

c dl 


(17) 


where  V = {Fv,  F);  ,o}  is  the  velocity  of  the  SCM  mass  centre;  0)  = {0,  0,  0)}  is  the  angular  velocity; 
m is  the  mass  of  model;  Ic  is  the  moment  of  inertia  of  model  about  mass  centre. 
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Passing  on  to  the  integration  along  the  longitudinal  absolute  co-ordinate  x and  adding  obvious 
kinematic  relations,  we  obtain  the  calculation  set  of  five  ordinary  differential  equations 


V cos (y/-a)^  = coV  +—1FX, 


dx 


m 


(18) 


, x d V v 1 

V cos(y/-cc) — — = -coVx  + — EF  , 
dx 


m 


(19) 


V cos(y/-a)—  = — ZM„, 
dx  L 


(20) 


V cos(t j/-a)^-  = (0, 
dx 


(21) 


— = tan(y/-a). 
dx 


(22) 


Here,  y/  is  the  model  pitch;  a is  the  angle  of  attack;  y is  the  ordinate  of  the  SCM  mass  centre  in  the 
absolute  co-ordinate  system  (see  Fig.  17). 


4.4.  Relations  for  forces  and  moments 


Two  kinds  of  forces  in  right  parts  of  Eqs.  (18)  - (20)  are  essential  in  considered  velocity  range  (see 
Fig.  17): 

- hydrodynamic  force  and  moment  on  the  cavitator  Fn , M n ; 

-hydrodynamic  force  and  moment  due  to  interaction  between  the  SCM  tail  and  internal  cavity  wall 

FX,MX. 

Coefficients  of  the  force  projections  acting  on  the  cavitator-disk  inclined  to  the  stream  at  angle  a are 
approximately  calculated  by  formulae  [2]: 

cx  - cx0cos2  a,  = cv0sina  cos2  a.  (23) 


For  disk  cavitators  with  plane  perpendicular  to  the  model  axis,  the  vector  Fn  is  always  directed  along 
the  model  axis.  Therefore,  the  force  on  the  cavitator  does  not  create  any  moment.  Formulae  (23)  agree 
with  experiment  for  a disk  and  blunted  cavitators  when  a < 50° . 

We  use  the  approximate  method  [9]  to  calculate  the  hydrodynamic  force  on  inclined  non-disk 
cavitator.  It  is  based  on  the  exact  solution  of  the  two-dimensional  problem  on  free  jet  flow  around  an 
inclined  polygonal  contour. 
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Process  of  interaction  between  the  SCM  tail  and  internal  cavity  wall  is  considered  as  unsteady 
planing  of  a prolate  body  along  the  curvilinear  liquid  boundary.  The  forces  arising  in  this  case  are 
calculated  basing  on  the  hypothesis  of  plane  sections  [2],  The  solution  of  the  problem  on  immersion  of 
the  circular  arc  through  the  curvilinear  boundary  is  used  in  this  case  [20]. 

If  both  the  middle  clearance  A = Rc  - Rs  and  the  model  immersion  depth  into  the  cavity  surface  hk 

are  small,  formula  for  the  cross  component  of  the  tail  force  in  the  body  coordinate  system  was  obtained 
in  the  form  [20] 


Fsy  = pnRs  V 


V , 


h(2  + h)  _2h_ 

(l  + hf  2j  + h_’ 


(24) 


where  Rs  is  the  stern  section  (transom)  radius  of  the  model;  l'\  is  the  cross  velocity  of  the  model 
transom;  V2  = -3  Rc  / dt  is  the  velocity  of  the  cavity  boundary; 
h = -hk  / A,  k = 1.2;  h\  < 0,  h2  < 0 is  the  model  transom  immersion  into  the  lower  and  upper  cavity 
boundaries,  respectively. 

The  longitudinal  component  of  force  Fsx  has  viscous  nature  and  is  calculated  by  formulae 


Fsx=^—Swcf{  Re), 


(25) 


where  Sw  is  the  area  of  washed  part  of  the  model;  c / is  the  viscous  drag  coefficient;  Re  is  the 
Reynolds  number. 

If  the  motion  velocity  is  very  high  V > 1000  m/s,  then  aerodynamic  forces  of  interaction  between  the 
model  body  and  vapour-gas-spray  medium  filling  the  cavity  can  essentially  influence  on  the  SCM 
dynamics  [1].  Our  evaluations  have  shown  that  this  influence  has  damping  character  [20]. 


4.5.  Stability  of  SCM  motion 


Computer  simulation  with  the  STAB  software  allows  to  investigate  the  stability  of  SCM  motion  "as 
a whole",  when  the  model  oscillates  within  the  cavity  ricocheting  by  its  tail  from  lower  and  upper  cavity 
walls  in  turn. 

The  PC-screen  view  after  execution  of  the  function  “Motion”  is  shown  in  Fig.  18,  a.  Here,  HX  is  the 
dimensionless  calculation  step;  N = 10  is  the  number  of  the  model  contacts  with  upper  and  lower  cavity 
contours.  The  velocities  Vx,  V and  (O  are  output  in  dimensionless  by  L,  V0  form  (where  L is  length  of 
the  model).  The  angles  y/,  a are  output  in  radians.  A signs  of  the  coefficients  show  the  direction  of  the 
force  and  moment  action. 

The  PC-screen  view  after  execution  of  the  function  “History”  for  the  same  model  and  starting  data  is 
shown  in  Fig.  18,  b.  The  y-coordinates  of  the  graphs  are  plotted  in  dimensionless  form,  the  x-coordinate 
x is  given  in  meters.  As  it  is  visible  from  graphs,  the  interaction  of  the  SCM  tail  and  water  has  behaviour 
of  short-term  impulses.  It  is  displayed  in  discontinuous  behaviour  of  functions  CO  (x)  and  V y (x)  [16]. 
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The  performed  analysis  has  shown  that  the  motion  stability  of  the  SCM  with  a disk  cavitator  is 
determined  by  three  dimensionless  parameters: 


/ 


c 


St  = 


COqL 

Ro  ’ 


A = 


Rc-Rs 

Rn 


Value  of  the  parameter  A is  determined  mainly  by  cavitator  radius  Rn  for  given  model  shape  and 
weakly  depends  on  starting  velocity  V0  . 

After  a number  of  calculations  at  constant  values  of  F0 , Ic , x2  and  varied  values  of  Rn , Q)q  , the 
plane  of  the  parameters  A,  St  is  divided  on  the  zone  I where  motion  is  stable  and  the  zone  II  where  the 
motion  is  unstable.  For  dynamically  similar  models  Ic  = const  the  stability  zones  are  equal. 

The  motion  stability  loss  is  visually  perceived  on  the  PC-screen  as  washing  the  frontal  part  of 
model  and  inadmissible  great  increase  of  the  pitch  y/. 

Also,  the  STAB  software  permits  to  simulate  the  regime  of  continuous  planing  the  model  along  the 
lower  cavity  wall  including  the  ventilated  one  and  with  working  propulsor.  In  this  case,  the  cavitator  is 
fixed  at  the  angle  of  8 to  the  model  axis  to  compensate  the  moment  of  the  force  arising  on  the  tail. 
Thrust  of  the  propulsor  is  added  to  the  right  part  of  Eq.  ( 1 8). 

The  computer  simulation  that  the  SCM  motion  in  regime  of  continuous  planing  is  unstable  when  the 
angle  8 is  constant. 


4.6.  Comparison  with  experiment 

We  confirmed  reality  of  the  SCM  self-stabilisation  mechanism  by  ricocheting  from  the  cavity  walls 
by  direct  shooting  the  motion  in  water  with  velocities  up  to  1200  m/s  [4,  6].  In  experiment,  its  action 
appears  in  periodic  supercavity  surface  perturbations  which  develop  next  according  to  the 
"independence  principle".  In  Fig.  19,  the  experimental  photographs  of  the  model  in  moment  of  touch  to 
the  cavity  wall  (compare  with  Fig.  14,  a)  and  the  supercavity  part  perturbed  after  contact  with  model  are 

shown.  In  this  case,  the  exposure  time  of  shooting  was  3 • 1 0 fl . 

Asymmetric  washing  the  model  occurs  when  the  SCM  motion  loses  stability.  It  results  in  impact 
increase  of  the  hydrodynamic  drag  and  deformation  of  model. 


5.  Schemes  of  the  hydrodynamic  stabilization  and  control  of  SC  objects 


5.1.  Scheme  with  free  planing  along  a cavity 


This  scheme  (Fig.  20)  of  the  vehicle  motion  is  more  favorable,  because  has  the  smallest  washed  part 
of  the  body  and,  hence,  the  lowest  hydrodynamic  drag.  However,  this  scheme  of  motion  will  have 
essential  distinctions  in  creation  of  stabilizing  hydrodynamic  forces. 
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At  motion  the  vehicle  is  located  in  the  cavity  under  some  angle  of  attack aA.  It  stipulates  the 
appearance  of  lift  on  the  tail  part  due  to  asymmetric  washing  the  cylindrical  part  of  the  body.  At  motion 
with  constant  velocity  and  depth  H (or  the  pressure  P0  = pgH ) we  can  choose  the  constructive 
parameters  of  the  vehicle: 

- position  of  gravity  center; 

- cavitator  dimensions; 

- deflection  angle  of aP; 

- gas-supply  value  of  the  artificial  cavity  Q or  its  length  Lc  = F(Q)  so  that  the  basic  conditions  of 

the  motion  stability  (equality  of  projection  of  all  forces  and  moments  acting  on  the  body  in  the  uniform 
rectilinear  motion)  are  fulfilled: 


Yc  + YB  =G-T  sin  a A , 

Xc  + X B - T cos«  4 , (26) 

M Q + Mg  = 0. 


The  hydrodynamic  forces  on  the  cavitator  at  deflection  on  small  angle  cxc  can  be  determined  by 
formulae  [13] 


YC  = XC0  cosac  sin  ac, 
Xc  = Xco  cos2  ac, 


(27) 


where  XCq  is  the  cavitator  drag  at  ac  = 0.  For  a disk: 


pV 

Xco=0.82^—S0, 


where  iSgis  the  area  of  the  cavitator  section  on  the  cavitating  edge;  OCc  =(Xp  —C/.A  is  the  cavitator 
angle  of  attack;  aP  is  the  positional  angle  of  the  cavitator  fixation;  CXA  is  the  body's  angle  of  attack; 
A B is  the  drag  of  washed  part  of  the  body  which  includes  the  friction  drag  of  the  side  surface  X F . 

At  this  flow  scheme  it  is  possible  to  save  the  traditional  principle  of  stabilization  and  control  due  to 
the  tail  fins  and  rudders,  but  with  an  additional  condition  - the  cross  stability  preservation. 

Really,  at  presence  of  roll  angle,  y,  the  lift  on  the  cavitator  deflected  on  angle  ac  will  form  the 
side  force  Zc  = Yc  sin  y , which  causes  to  the  body  yaw  (Fig.  21). 

The  cross  stabilization  of  the  object  can  be  obtained  by  following  ways: 
a)  due  to  the  location  of  mass  center  lower  than  a line  connecting  the  points  Yb  of  application  of 
the  forces  (Fig.  22).  The  stabilizing  moment  My  arises  at  arising  the  roll  y,  (Fig.  22).  In  this  case  it  is 
possible  to  stabilize  the  object  without  fins  at  all  at  sufficient  magnitude  of  My 


My  = Ghc  sin  y ; 


(28) 
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b)  due  to  the  hydrodynamic  fins  of  roll  (Fig.  23)  creating  the  stabilizing  moment. 

A possible  scheme  of  cross  stabilization  of  the  body  within  the  cavity  due  to  two  fins  is  shown  in 
Fig.  23. 

There  the  tail  part  of  the  vehicle  is  equipped  by  two  supercavitating  fins  creating  hydrodynamic 
lifts:  the  right,  and  left,  Y2  . At  possible  angle  of  roll  of  the  vehicle  on  angle  \gamma  the  stabilizing 
moment  My  must  arise  due  to  a difference  of  these  forces  caused  by  various  immersion  of  the  right 
and  left  fins  into  flow  through  the  free  cavity  boundary.  The  total  lift  of  the  two  fins  must  not 
exceed  10  4-15  of  the  total  lift,  YB , on  the  vehicle  tail 


r1+r2=o.i57g. 

The  hydrodynamic  forces  on  fins  in  stable  position  arise  under  action  of  the  hull  slope  about  the 
cavity  axis  at  gravity  influence  and  also  flow  bend  in  the  back  cavity  part  caused  by  cavity  floating-up 
effect. 


Stabilization  in  horizontal  plane  is  achieved  due  to  using  a number  of  factors: 

- stabilizing  influence  of  the  hydrodynamic  forces  of  washed  part  of  the  hull  located  past  the  mass 
center  of  supercavitating  object  (SCO); 

- application  of  neutral  stability  cavitators  (disk)  relatively  to  angles  of  deviation  in  the  horizontal 
plane; 

- using  the  schemes  of  a self-stabilization  by  roll: 

- low  location  of  the  mass  center, 

- V-shaped  system  of  fins. 


Stabilization  of  SCO  in  vertical  plane  is  reached  due  to  the  following  factors: 

- the  cavitator  with  deviation  angle  having  stability  about  angular  deviations  in  vertical  plane  (disk,  cone 
P > 110° ) [14]; 

- stabilizing  influence  of  the  washed  hull  part  located  past  the  SCO  mass  center; 

- stabilizing  influence  of  the  V-shaped  system  consisting  of  two  fins  located  on  the  tail  end  of  washed 
part  of  SCO. 


The  vehicle  completely  placed  in  the  cavity  is  shown  in  the  scheme  in  Fig.  20.  The  supporting 
hydrodynamic  forces  are  created  by  deflecting  statically  stable  cavitator  [14],  Yc,  and  planing  tail  part 
of  the  vehicle,  YB . 

The  low  location  of  a point  of  the  hydrodynamic  force,  Yb>  application  hampers  the  execution  of 
cross  self-stabilization  scheme  by  lowing  the  mass  center  me  below  than  a line  connecting  the 
points  of  the  hydrodynamic  force  application  on  the  cavitator,  Oc,  and  on  the  tail,  0B , of  the  body  (it  is 
shown  by  dotted  line  in  the  scheme). 

The  considered  stabilization  schemes  do  not  exclude  the  application  of  means  of  the  active  control 
of  vehicle  for  it  stabilization  during  motion. 
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For  SCO  control  may  be  used  general  aerohydrodynamics  principle  - to  realize  the  object  control  by 
using  elements  creating  hydrodynamic  side  forces  as  wings,  flaps,  rudders,  jets. 

In  SC  flow  scheme  the  cavitator  is  a convenient  element  for  side  force  creation  owing  formulae  (27). 

SCO  control  in  horizontal  plane 

• by  cavitator  deflection  in  horizontal  plane; 

• by  cavitator  rotation  in  corresponding  side  according  scheme  (Fig.  21); 

• by  SCO  hull  rotation  in  corresponding  side  by  using  active  fins  (Fig.  23). 

SCO  control  in  vertical  plane 

• by  cavitator  deflection  in  vertical  plane  according  scheme  in  Fig.  20  and  formula  (27). 


5.2.  Optimization  of  fin  number 

Traditional  stabilization  schemes  of  bodies  of  revolution  with  large  aspect  ratio  require  presence  of 
two  pairs  of  fins: 

- the  first  pair  to  stabilize  motion  in  horizontal  plane; 

- the  second  pair  to  stabilize  motion  in  vertical  plane  [4], 

Aerodynamic  schemes  for  stabilization  and  control  were  developed  to  minimize  the  drag  of  aircraft 
and  rockets.  They  have  three  elements  located  by  120°  or  T -shaped  and  also  two  elements  at  V - 
shaped  location  [4]. 

The  most  of  the  hydrodynamic  high-speed  vehicles  of  known  designs  have  two  pairs  of  fins  and 
rudders  for  the  motion  stabilization  reliability,  control  convenience  [6]  including  supercavitating  ones 

[9]- 


5.3.  Optimization  of  fin  shape 

Use  of  the  rudders  and  fins  on  supercavitating  objects  has  essential  features: 

1)  two-phase  regime  of  flow  around  the  part  or  whole  fin  installed  past  the  artificial  supercavity; 

2)  incomplete  washing  the  fin  at  possible  air  layer  on  the  supercavitating  object  hull  in  the  fin  mounting 
place; 

3)  high  overloads  at  water  entry  and  take  off  of  the  fins  at  oscillations  of  the  cavity  boundary  or  the 
body  in  the  supercavity. 

Presence  of  the  two-phase  flow  promotes  the  flow  separation  from  the  frontal  edge  of  the  fin  and 
cavity  formation.  It  results  in  abrupt  reduction  of  the  fin  lift. 

The  air  layer  presence  on  the  hull  in  the  fin  mounting  place  also  causes  to  the  lift  loss  not  only  due 
to  the  loss  of  the  fin  interference  with  the  hull.  In  the  case  of  the  plane  wall  this  loss  is  equal  to 

reduction  of  the  effective  aspect  ratio  A in  two  times  according  to  formula  for  : 


c° 


2k 


(29) 
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High  overloads  at  abrupt  change  of  medium:  water-gas-water  require  increased  strength  and  increasing 
the  fin  profile  thickness. 

All  these  features  of  flow  around  the  rudders  and  fins  on  the  supercavitating  body  of  revolution  were 
investigated  theoretically  and  experimentally  at  the  IHM  UNAS.  These  research  results  were 
development  of  two-phase  profile  stable  for  variable  loads  [12,13]  and  recommendations  for  the  fin 
shape  [14]. 

All  the  rudders  has  wedge  profile  of  the  cross  section  {fi  = 3°  -4-  6 ) (Figs.  24,  25). 

The  experimental  results  show  that  the  wedge  profile  gives  more  stable  flow  characteristics  for 
regimes  of  presence  of  gas  and  water-air  layers  on  the  body  hull.  In  this  case  the  fins  of  small  aspect 
ratio  A < 1 are  more  preferable,  because  they  have  no  separations  in  great  range  of  angles  of  attack. 
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